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A  numerical  model  was  implemented  to  analyze  the  thermodynamic  performance  of  the  co-  and  counter¬ 
flow  operations  of  an  anode-supported  direct  internal  reforming  (DIR)  planar  solid  oxide  fuel  cell  (SOFC). 
This  developed  model  was  validated  by  comparing  with  experimental  and  simulated  results  taken  from 
the  literature.  The  model  is  capable  of  capturing  the  detailed  distribution  of  the  local  temperatures,  species 
concentrations,  current  density,  and  polarization  losses  in  streamwise  direction.  Energy  and  exergy  con¬ 
cepts  were  used  to  evaluate  the  DIR-SOFC  performance  under  co-  and  counter-flow  operations.  The  study 
indicates  the  energy  and  exergy  efficiencies  of  DIR-SOFC  performance  under  co-flow  operation  are  more 
sensitive  to  the  increase  of  current  density  than  that  under  counter-flow  operation.  Particular  attention 
was  paid  to  cell  temperature  profiles  to  avoid  mechanical  failure  due  to  high  thermal  stresses.  The  result 
shows  that  the  material  constraints  need  to  be  considered  as  well  as  the  energy  and  rational  efficiencies 
in  evaluating  the  performance  of  SOFC.  The  preferred  flow  configuration  can  be  changed  depending  on 
the  cell  geometry  and  operation  conditions  if  we  consider  the  material  constraints. 

©  201 1  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Among  fuel  cells,  high  operating  temperature  SOFCs  could 
improve  and  promote  the  exploitation  of  power  generation  system 
as  combined  heat  and  power  systems  to  raise  the  energetic  effi¬ 
ciency  [1-7].  Furthermore,  high  operating  temperature  SOFCs  also 
offer  many  advantages  over  other  fuel  cells  (i.e.,  fuel  flexibility,  tol¬ 
erance  to  impurities).  However,  this  high  operating  temperature 
poses  a  problem  in  terms  of  the  thermal  management  of  SOFCs  [8]. 
Because  SOFCs  are  made  of  fragile  ceramic  materials,  a  large  tem¬ 
perature  gradient  may  cause  cracking  of  the  cell.  It  is  required  to 
properly  control  the  cell  temperature  distribution  for  its  safe  and 
efficient  operation. 

One  of  the  attracting  features  of  SOFC  is  the  direct  internal 
reforming  (DIR)  of  hydrocarbon  fuels.  The  endothermic  reforming 
reaction  proceeds  on  the  Ni  catalyst  in  the  anode  in  the  direct  inter¬ 
nal  reforming  process.  It  utilizes  the  waste  heat  generated  by  the 
electrochemical  reaction  and  other  irreversible  processes  to  offset 
the  heat  requirements  of  the  reforming  reaction,  resulting  in  the 
increase  in  the  performance  of  the  SOFC.  It  also  reduces  the  amount 
of  the  extra  air  supplied  to  the  cathode  channel  because  the  cool¬ 
ing  demand  is  reduced,  thus  lowering  the  cost.  One  drawback  is 
that  it  may  cause  a  large  temperature  gradient  to  the  cell  because 
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of  the  strong  endothermic  nature  of  the  reforming  reaction.  How¬ 
ever,  direct  measurements  of  the  temperature  distribution  in  an 
operating  DIR-SOFC  are  difficult  because  of  the  high  temperature, 
narrow  channel  gaps,  and  small  flow  rates. 

Numerical  simulation  is  expected  to  be  an  effective  tool  to 
understand  the  phenomena  in  a  DIR-SOFC  system.  CFD-based 
two-  and  three-dimensional  simulations  of  fuel  cells  can  be  found 
in  the  literature  [9-11]  but  owing  to  the  heavy  computational 
load,  the  computational  domain  is  generally  limited  to  a  single 
module  or  channel.  Iwai  et  al.  [12]  recently  developed  a  CFD 
based  quasi-three-dimensional  simulation  model  that  applies  the 
volume-averaging  method  to  the  flow  passages  by  assuming  that 
a  porous  material  is  inserted  as  a  current  collector.  It  reduces  the 
computational  time  and  cost  while  maintaining  the  ability  to  solve 
the  flow  and  pressure  fields  in  the  SOFC.  The  CFD  based  simulation 
is  effective  to  investigate  detailed  phenomena  in  the  SOFC  but  it 
is  an  expensive  computational  process  for  a  system  level  analysis 
where  modules  are  usually  described  with  lumped  (O-dimensional) 
models.  The  use  of  the  lumped  model  to  DIR-SOFC,  however,  needs 
to  be  carefully  examined  before  its  application  to  a  system  analy¬ 
sis  because  the  large  temperature  gradient  expected  in  a  DIR-SOFC 
and  its  effects  cannot  be  taken  into  account  if  a  lumped  model  is 
used.  To  consider  the  distributions  of  variables  in  the  DIR-SOFC,  1  -D 
model  is  at  least  needed. 

Campanari  [13]  proposed  two-dimensional  CFD  based  numer¬ 
ical  models  for  electrolyte-supported  DIR-SOFC.  They  performed 
parametric  analysis  on  the  effects  of  heat  losses,  air  ratio  and  cell 
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Nomenclature 

A  cell  active  area,  m2 

Ac  cross  section  area,  m2 

As  surface  area,  m2 

Cp  specific  heat  at  constant  pressure,  J  kg-1  K 

Di  eff  effective  diffusivities  of  the  species  z,  m2  s-1 

Dim  diffusivity  of  the  species  z,  m2  s-1 

Dik  Knudsen  diffusion  of  the  species  z,  m2  s-1 

dh  channel  hydraulic  diameter,  m 

dpore  pore  diameter,  p,m 

ex  specific  exergy  flow,  J  kg-1 

Ech4  activation  energy  of  the  methane  reforming  reaction 

Ex  exergy  flow  rate,  W 

F  faraday  constant,  96487  C  mol-1 

h  convective  heat  transfer  coefficient,  W  m-2  K 

ic  current  density,  A  nrr 2 

i0  exchange  current  density,  A  nrr2 

kshift  rate  constant  of  forward  shift  reaction, 

mol  nr 3  Pa2  s 
I<  permeability,  m2 

Kp  equilibrium  constants 

zrz  mass  flow,  kgs-1 

Mi  molecular  weight  of  species  i,  kg  mol-1 
ne  number  of  electrons  participating  in  the  electro¬ 

chemical  reaction 

rij  mole  flow  rate  of  species  z,  mol  s-1 

P  pressure,  Pa 

Peiec  electrical  power,  W 

Pi  partial  pressure  of  species  z,  Pa 

Q.  heat  generation  rate,  W  nrr3 

R reform  reforming  reaction  rate,  mol  nr3  s 
Rshift  shift  reaction  rate,  mol  nrr3  s 

R  universal  gas  constant,  8.314  J  mol-1  K 

T  temperature,  K 

V  electric  potential,  V 

V0c  open-circuit  voltage,  V 

Greek  letters 

A G°  change  of  standard  Gibbs  free  energy,  J  mol-1 

AH  enthalpy  change,  J  mol-1 

sp  porosity 

rj  energetic  efficiency 

Pact  acticvation  polarization  loss,  V 

Pconc  concentration  polarization  loss,  V 

Pohm  ohmic  polarization  loss,  V 

X  thermal  conductivity,  W  nrr1  K 

Pair  air  stream  density,  kg  nr 3 

Pfuei  fuel  stream  density,  kg  nrr3 

Pj  specific  electrical  resistance  of  the  component  j, 

C2m2 

< j  ionic/electronic  conductivity,  S  nr1 

8  thickness,  m 

t  tortuosity 

x/f  rational  efficiency 

Subscripts 

0  properties  of  the  environment 

air  air,  air  channel 

an  anode 

ca  cathode 

ch  chemical 

co-flow  co-flow  cell  configuration 

counter-flow  counter-flow  cell  configuration 


fuel  gas  mixture  at  the  fuel  channel,  fuel  channel 

H2  hydrogen  oxidation 

in  inlet  value 

out  outlet  value 

PEN  positive-electrolyte-negative  structure 

ph  physical 

shift  shift  reaction 

SOFC  fuel  cell 

SR  steam  reforming 

TPB  three-phase  boundary. 


configurations  (cross-,  co-  and  counter-flow  cells)  on  the  tem¬ 
perature  distributions  of  a  planar  electrolyte-supported  DIR-SOFC. 
To  improve  power  output,  an  anode-supported  planer  DIR-SOFC 
recently  gains  more  attentions.  Aguiar  et  al.  [14]  and  Li  et  al. 
[15]  developed  one-dimensional  model  of  an  anode-supported 
planer  DIR-SOFC  and  applied  it  to  the  co-  and  counter-flow 
operations.  Both  reported  the  detailed  distributions  of  the  cell 
temperature,  species  more  fractions,  current  density  and  elec¬ 
trochemical  performance.  Flo  et  al.  [16]  performed  CFD-based 
three-dimensional  simulations  using  a  commercial  package  (Star- 
CD)  and  reported  similar  results.  All  these  studies  [14-16]  agree 
that  the  counter-flow  configuration  has  advantages  with  high  per¬ 
formance  compared  to  the  co-flow  configuration.  However,  an 
unfavourable  steep  temperature  gradient  is  also  commonly  pre¬ 
dicted  for  the  counter-flow  configuration.  Ho  et  al.  reported  that  a 
decrease  of  the  air  inlet  temperature  alleviated  the  large  tempera¬ 
ture  gradients  in  the  counter-flow  case  but  it  inevitably  decreases 
the  cell  average  temperature  and  thus  lowers  the  cell  efficiency. 
In  a  DIR-SOFC  where  a  non-uniform  temperature  distribution  is 
expected,  the  selection  of  the  configurations  and  operation  condi¬ 
tions  need  to  be  done  considering  not  only  the  efficiency  but  the 
material  constraints. 

In  this  paper,  we  develop  a  numerical  model  of  an  anode- 
supported  planer  DIR-SOFC  to  clarify  suitable  operating  parameters 
for  co-  and  counter-flow  operations  considering  both  the  effi¬ 
ciency  and  material  constraints.  With  its  future  application  to  a 
system  analysis  in  mind,  the  model  is  one-dimensional  and  capa¬ 
ble  of  performing  exergy  analysis.  A  comparative  performance 
study  between  co-  and  counter-flow  planar  anode-supported  DIR- 
SOFCs  is  performed  under  constant  fuel  and  air  utilization  factors 
with  the  distribution  of  the  temperatures,  species  concentrations, 
current  density,  and  polarization  losses.  Since  changing  inlet  tem¬ 
peratures  can  have  either  a  beneficial  or  an  undesirable  impact  on 
performance  and  life  span  of  the  cell,  effects  of  the  current  density 
and  inlet  temperatures  under  co-  and  counter-flow  operations  are 
investigated  by  energy  and  exergy  methods  to  clarify  the  maximum 
cell  performance  with  lowest  risk  of  thermal  failure.  The  previous 
studies  in  this  category  found  in  literature  report  the  advantage 
of  a  counter-flow  configuration  over  a  co-flow  configuration  from 
the  energy  efficiency  point  of  view  [14-16].  In  this  study  we  revisit 
the  discussion  on  the  preferred  flow  configuration  with  consider¬ 
ation  of  the  material  constraints  as  well  as  the  energy  and  exergy 
efficiencies. 

2.  Model  description 

2.1.  Computational  domain 

In  a  practical  planer  SOFC  system,  a  number  of  cells  are  stacked 
in  series.  When  the  number  of  the  stacked  cells  is  sufficiently  large, 
most  of  the  cells  located  in  the  core  of  the  stack  are  similarly  sur¬ 
rounded  by  other  cells  and  are  therefore  operated  in  a  similar 
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situation  to  each  other.  Attention  is  therefore  focused  in  this  study 
to  a  single  cell  in  such  a  situation  in  a  stack  assuming  that  it  repre¬ 
sents  the  entire  stack  characteristics.  In  other  words,  we  assume  an 
ideal  thermal  insulation  of  the  system  and  the  cells  in  the  stack  are 
identical.  As  shown  in  Fig.  1,  a  planar  SOFC  of  co-  and  counter-flow 
types  was  modeled  with  the  pair  of  two  channels  under  a  steady- 
state  condition;  the  SOFC  is  composed  of  the  fuel  and  air  channels, 
the  PEN  structure  (anode,  electrolyte  and  cathode)  and  intercon¬ 
nection.  As  shown  in  Fig.  2,  heat  loss  near  the  edge  of  the  cell  stack 
is  also  neglected. 

The  electrochemical  reaction  is  considered  to  be  attributed  to 
only  hydrogen;  the  electrochemical  fuel  value  of  CO  is  readily 
exchanged  for  hydrogen  by  the  rapid  shift  reaction  assuming  chem¬ 
ical  equilibrium  [17].  In  other  words,  CO  is  considered  to  take  part 
only  in  the  shift  reaction.  Whereas  the  open-circuit  voltage  ( V0c)  of 
an  internal  reforming  SOFC  is  calculated  according  to  the  electro¬ 
chemical  oxidation  of  H2,  the  species’  consumption  and  production 
is  determined  collectively  from  the  steam  reforming  reaction  and 
shift  reaction. 

As  presented  in  Fig.  2,  the  one  channel  region  is  discretized 
into  90  control  volumes  of  uniform  size  (AX)  in  the  x-direction. 
Each  control  volume  contains  three  temperatures  corresponding 
to  air  and  fuel  streams,  Tj fuel  and  Tair,  and  temperature  of  the  cell 
or  the  solid  part  (PEN  and  interconnect  are  lumped  together),  Tso/ld, 
respectively.  The  mass  and  heat  balances  of  each  control  volume  are 
analyzed  and  solved  by  finite  difference  method.  Fig.  2  also  shows 
that  the  fuel  always  flows  in  the  x-direction  while  air  flows  either 
in  the  x-direction  (co-flow  operation)  or  in  the  negative  x-direction 
(counter-flow  operation). 


2.2.  Mass  balances 


Mass  balances  are  formulated  for  each  species  on  the  basis 
of  the  relationship  between  the  local  current  and  the  change  in 
the  concentrations.  The  following  three  equations  summarizes  the 
reactions  considered  in  the  cell. 


Steam  reforming  reaction  :  CH4  +  H20  CO  +  3H2  ( 1 ) 

Shift  reaction  :  CO  +  H20  C02  +  H2  (2) 

Overall  cell  reaction  :  FI2  +  ^02  FI20  (3) 

The  open-circuit  voltage  ( V0c).  is  described  by  the  Nernst  equa¬ 
tion  as  a  function  of  operating  temperature  (T)  and  partial  pressure 

(p). 


Voc  = 


A  G°  RT 
IT  +  2 F 


•  In 


Ph2  ■Pof/2\ 
Ph2o  ) 


(4) 


These  reactions  lead  to  five  mass  balances  over  the  unit  cell 
where  x,  y  and  z  are  the  conversion  rates  of  CFI4,  CO  and  H2,  respec¬ 
tively. 


AnCH4  =  -x 

(5a) 

Anco  =x-y 

(5b) 

AnCo2  =  y 

(5c) 

Ahh2  =3  x  +  y-z 

(5d) 

Ahh2o  =  -x-y  +  z 

(5f) 

Air  flow 
(coflow) 

Air  flow  4 - 

(counter  flow) 


ar 


Sc 


=  0 


*-0 

Fuel  flow 


T solid 


1  fuel 


Ax 


interconnect 


air  channel  „ 


PEN 


fuel  channel 


interconnect 


Air  flow 
(co-flow) 

— ► 


S71 


Sc 

— ► 


=  0 


x=L 


Fig.  2.  Discretized  domain  for  co-  and  counter-flow  configuration. 
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A  model  proposed  by  Achenbach  [  1 8  ]  is  adopted  for  the  reaction 
rate  of  the  methane  reforming  reaction: 


rCH4  =  kcH4PCH4 


(6) 


where,  Ech4  (=  82  kj  mol-1 )  is  the  activation  energy  of  the  reaction, 
/<ch4  (=  4274  mol  m_2-bar-s)  is  the  pre-exponential  factor  and  pCH4 
is  the  partial  pressure  of  methane. 

The  chemical  equilibrium  constant  of  the  shift  reaction  is 
another  constraint  governing  relations  among  chemical  compo¬ 
nents  in  the  control  volume. 


Kp, shift 


Pc o2  -Ph2 
Pc o  •  Ph2o 


(nc  o2  +y/ntQt  +  2x)(nH2  +  3x  +  y  -  z/ntot  +  2x) 
(nco  +  x  -  y/ntot  +  2x)(nH2o  -x-y  +  z/ntot  +  2x) 


(7) 


Here,  rp  denotes  the  molar  flow  rate  of  species  i  at  the  inlet  of 
the  control  volume.  ntot  is  the  total  molar  flow  rate  of  all  chem¬ 
ical  components  at  the  same  location.  If  the  exit  temperature  of 
each  control  volume  is  evaluated  from  energy  balance  equations 
described  in  Section  2.4,  chemical  equilibrium  constant  (7)  can  be 
obtained  and  applied  to  analyze  mass  conservations  of  chemical 
components  in  each  control  volume.  As  mentioned  above,  the  shift 
reaction  is  always  assumed  to  be  at  equilibrium  in  this  study.  The 
equilibrium  constant  can  be  given  by  (8)  [19]: 

I<pshift  =  exp(-0.2935 Z3  +  0.6351 Z2  +  4.1788Z  +  0.3619)  (8) 

where  Z=  ( 1 000/T)  - 1 . 


Table  1 

Calculation  based  parameters  for  diffusion  overvoltage. 


Parameter 

Value 

Electrode  porosity  ( ep ) 

0.48 

Electrode  tortuosity  (r) 

5.4 

Pore  diameter  (dpore)  (p,m) 

1.0 

Taken  from  [22]. 

where  sp  and  r  are  the  porosity  and  tortuosity  of  electrode  mate¬ 
rials,  respectively,  shown  in  Table  1.  Diffusivity  of  species  i  in 
multicomponent  gas  mixture,  Dim  can  be  estimated  by  Wilke  [21  ]: 


Dim  = 


1  -Xj 


±  Mj/Dij ) 


(15) 


where  X;  is  the  mole  fraction  of  gas  species  i.  For  prediction  of  the 
binary  diffusivity  of  the  gas  mixture  composed  of  species  i  and  j 
(Dy ),  the  Fuller-Schettler-Giddings  formula  [20]  is  adopted  in  this 
study: 


Dy  = 


0.1013r175((l  /M, )  +  (1  /JW2))° 


(E"i)1/3  +  (E"2 


4/3' 


(16) 


where  is  molecular  weight  of  species  i;  (Eu*)  represents  diffu¬ 
sion  volume  of  species  i,  which  values  can  be  found  in  [20],  Knudsen 
diffusion  ( Dik )  can  be  expressed  as: 

/  T  \0,5 

Dik  =  48.5dPore(wJ  (17) 


2.3.  Electrochemical  model 


where  dpore  is  pore  diameter. 


2.3.1.  Concentration  polarization 

The  concentration  polarization  is  ascribed  to  the  finite  gas  dif¬ 
fusion  processes  that  govern  movement  of  gasses  into  and  out  of 
the  electrochemical  reaction  site.  The  potential  difference  between 
operations  with  and  without  current  is  the  concentration  polariza¬ 
tion,  hence  is  evaluated  as: 


P  cone,  an 

ceb 

ii 

H 

(  Ph2o,tpbPh2  \ 

^Ph2oPh2,tpb  J 

(9) 

Pconc,ca 

RT 

ln  ( 

(  P°2 

(10) 

~  ~4F~ 

vPo2,tpb  J 

where  the  subscript  TPB  denotes  the  three-phase  boundary.  To  cal¬ 
culate  the  pressure  at  the  reaction  sites,  the  following  equations 
have  been  adopted: 


Ph2,TPB  =  Pu2,an 


—  ic 


R  ■  TSan 

2 F  ■  DH2q  ejj 


Ph20,TPB  =  PH20,an  +  k 


R  ■  TSan 

2F  '  DH2,e ff 


(ID 

(12) 


P02,TPB  =  Pair  ~  ( Pair  ~  Po2,ca )  exp  (  ic  —  jj  ^ -  ]  (13) 

\  • 1  air  •  u02,eff  y 

In  the  equations,  P  represents  total  pressure,  and  Di  eff  represents 
the  effective  diffusivities  of  the  species  i  at  anode  and  cathode.  To 
evaluate  the  effective  diffusivities,  combined  ordinary  and  Knudsen 
diffusion  [20]  has  been  used: 


2.3.2.  Activation  polarization 

Chemical  reactions,  including  electrochemical  reactions, 
involve  energy  barriers  which  must  be  overcome  by  the  reacting 
species.  This  energy  barrier  is  called  the  ‘activation  energy’  and 
results  in  activation  or  charge-transfer  polarization,  which  is  due 
to  the  transfer  of  charges  between  the  electronic  and  the  ionic 
conductors  taking  place  at  the  three-phase  boundary  layer  of  both 
electrodes.  The  activation  polarization  may  be  regarded  as  the 
extra  potential  necessary  to  overcome  the  energy  barrier  of  the 
rate-determining  step  of  the  reaction  to  a  value  such  that  the  elec¬ 
trode  reaction  proceeds  at  a  desired  rate.  Activation  polarization 
(Pact)  is  normally  expressed  by  the  well  known  Butler  Volmer 
equation. 


ic  =  i0 1  exp  (  f  h'F^jCt^  ~  exp 


-(1-/1) 


PleFrjact  1  \ 

RT  J  J 


(18) 


The  parameters  ne  and  in  Eq.  (18)  were  set  equal  to  2  and  0.5, 
respectively  [14].  The  activation  overpotential  can  be  expressed  as 


Pact  — 


2  RT 
neF 


(19) 


Table  2 

Calculation  based  parameters  for  activation  overvoltage  in  Eqs.  (20)  and  (21). 


Parameter  Value 


Pre-exponential  factor  for  anode  ( yan )  ( A  m- 2 )  6.54  x  1 01 1 

Activation  energy  for  anode  (Eact,an )  0  mol-1 )  1 40,000 

Pre-exponential  factor  for  anode  ( ycat )  (Am-2)  2.35  x  1011 

Activation  energy  for  anode  ( Eactxat )  (J  mol-1 )  137,000 


Taken  from  [14,23]. 
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Table  3 

Specific  resistivity  for  the  model. 


Component 

Material 

Specific  resistivity  (pj)  (£2  m) 

Electrolyte 

YSZ 

2.94  x  10“5  exp(  10,350 /T)a 

Anode 

Ni/YSZ  cermet 

2.98  xlO-5  exp(-1392/T)a 

Cathode 

LSM-YSZ 

8.114  x  10“5  exp(600/T)a 

Interconnector 

Doped  LaCr03 

0.000321 5b 

a  Taken  from  [1]. 
b  Taken  from  [4]. 


where  i0  is  the  exchange  current  density.  Eqs.  (20)  and  (21)  are 
used  to  evaluate  the  values  of  the  exchange  current  density  for  the 
anode  and  the  cathode,  (see  variables  in  Table  2),  respectively. 

<—*(£)  -  (-%r)  <-> 

I).., %=)  (21) 

2.3.3.  Ohmic  loss 

The  ohmic  losses  in  SOFC  are  due  to:  (i)  the  electrons  flow 
through  the  anode,  cathode  and  interconnections;  (ii)  the  ionic 
flow  through  the  electrolyte.  The  resistivity  of  the  component  j  (pj) 
strongly  depends  on  the  temperature  [5]  as  presented  in  Table  3. 
Using  the  thickness  of  the  j-th  component,  Sj,  the  total  ohmic  loss 
(rj ohm)  can  be  expressed  as  follows: 

h ohm  =  Q  ^  (22) 

j 


Table  4 

Nusselt  number. 


Air  side 

Fuel  side 

Nu 

4.0 

6.2 

Taken  from  [28]. 


The  heat  generation  in  solid  structure  (Qsource, so  lid)  is  caused  by 
the  effect  of  electrochemical  reaction,  steam  methane  reforming 
reaction  and  the  ohmic  loss,  and  is  given  by  Eq.  (26),  where  AH  is 
the  change  in  enthalpy  of  each  reaction. 


Qsource, solid  —  *  x  AHsr  Z  X  Tso/lc/ ASh2  +  ic^hohm  (26) 

Heat  generation  associated  with  the  shift  reaction  occurs  in  the 
fuel  channel.  The  governing  energy  equation  for  the  fuel  channel 
is: 


0  —  PfuelCpjuel  AxAcjueiUfuei 


dTfuei 

dx 


hfuelAs,fuel(TSolid  7 fuel ) 


Qsource, fuel 


(27) 


The  heat  generation,  Qsource, fueh  in  the  fuel  channel  is  expressed 
as  follows: 

Qsource, fuel  —  x  shift  (28) 

2.5.  Exergy  related  to  a  stream  flow 


2.4.  Thermal  model 


The  streamwise  temperature  distributions  of  the  air  in  air  chan¬ 
nel,  Tair,  the  gas  mixture  in  fuel  channel,  ^fueb  and  solid  structure, 
T solid y  are  considered  in  the  present  1-D  numerical  simulation. 
The  heat  generation  and  heat  transfer  processes  considered  in  the 
model  are  the  heat  release  and  absorption  arising  from  the  electro¬ 
chemical  reaction,  the  electrical  resistance,  convective  heat  transfer 
between  the  solid  phase  and  gas  streams  and  the  conductive  heat 
transfer  in  solid  part.  At  an  elevated  temperature,  heat  transfer  by 
thermal  radiation  can  become  important.  However  it  is  reported  in 
literature  [24-27]  that  unlike  planar  electrolyze-supported  SOFC  or 
tubular  SOFC,  the  effect  of  radiation  heat  transfer  on  temperature 
profiles  is  insignificant  in  planar  anode-supported  SOFC  and  can  be 
negligible.  As  we  only  discuss  the  anode-supported  planer  SOFC  in 
this  study,  we  neglected  the  effects  of  thermal  radiation. 

Since  no  heat  generation  occurs  in  the  airside,  the  energy  con¬ 
servation  equation  in  air  channel  is  given  by: 


0  —  Pair^~P,air  A.xAcairUajr  ^  ^air^s,air(^solid  ^air)  (23) 


where  pair  is  the  density  of  gas  mixture  in  air  channel;  Cp  air  is  the 
specific  heat  capacity  of  gas  mixture  in  air  channel  and  Ac  air  is  the 
cross  section  area  of  the  channel.  The  heat  transfer  coefficient  h  is 
given  by  the  Nusselt  expression 


,  Nu  x  X 
h  =  —j — 

dh 


(24) 


where  <4  is  the  channel  hydraulic  diameter.  The  values  of  Nus¬ 
selt  number  and  thermal  conductivity  are  set  referring  literature 
[28,29]  and  are  summarized  in  Tables  4  and  5. 

For  the  solid  part  of  the  cell,  the  energy  conservation  equation 
is  given  by: 

q2j 

Qsource, solid  =  ^  solid  Ac,  solid  A.X  7^  +  ^ air As, air (7 solid  ~  4n'r) 


fuel  ^s,  fuel  ( T solid  Tfuel ) 


(25) 


Exergy  associated  with  material  stream  is  equal  to  the  maxi¬ 
mum  amount  of  work  obtainable  when  the  stream  is  brought  from 
its  initial  state  to  the  dead  state  by  processes.  Exergy  transfer  rate 
associated  with  material  stream  can  be  divided  into  physical  and 
chemical  exergy  components.  Physical  exergy  is  the  work  obtain¬ 
able  by  taking  the  substance  through  reversible  processes  from  its 
initial  state  temperature  T  and  pressure  P,  to  the  environmental 
state.  It  can  be  calculated  with 

E*ph  =  ^n,[(h  -  h0)  -T0(s-  s0)],  (29) 

i 

where  h,  represents  the  mole  flow  rate  of  species  i,  h  is  the  specific 
enthalpy  and  s  the  specific  entropy  and  the  properties  indicated 
with  the  subscript  0  refer  to  the  environmental  state. 

Chemical  exergy  is  equal  to  the  maximum  amount  of  work 
obtainable  when  the  substance  under  consideration  is  brought 
from  the  environmental  state,  defined  by  the  parameters  To  and 
P0,  to  the  reference  state  by  processes  involving  heat  transfer  and 
exchange  of  substances  only  with  the  environment.  The  chemical 
exergy  for  mixtures  can  be  calculated  as  follows: 

Exch  —  —  P'i,  00  )  (30) 

i 

where  /xl  0  is  the  chemical  potential  of  species  i. 

2.6.  Definitions  of  efficiencies 

In  energy  conversion  processes,  the  efficiencies  can  be  defined 
in  many  ways.  In  this  study  we  defined  two  kinds  of  efficiencies. 


Table  5 

Thermal  conductivity. 


Anode 

Cathode 

Electrolyte 

X  (Wm-1  K"1) 

11 

6 

2.7 

Taken  from  [29]. 
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Table  6 

Operating  parameters  and  SOFC  channel  geometries  used  in  electrochemical  model 
validation  [16]. 


Parameters 

Value 

Button  cell  geometries 

Anode  thickness  (pm) 

1000 

Cathode  thickness  (pm) 

20 

Electrolyte  thickness  (pm) 

8 

Fuel/air  stream  inlet  pressure  (bar) 

1 

Cell  mean  temperature  (I<) 

873-1073 

In  let  gas  composition 

Fuel: 

97%H2i  3%H20 

Air: 

21%  02  and  79%  N2 

2.6 A.  Energetic  efficiency  (rj) 

In  any  system,  energetic  efficiency  is  defined  as  the  ratio 
between  energy  in  product  outputs  to  the  energy  in  the  fuel  inputs. 
It  can  be  applied  to  SOFC  as  follows: 


i)  — 


_ Pelec _ 

Sch4,co,h 2(mfuel  x 


Fig.  3.  Comparison  of  the  simulation  results  with  experimental  data  by  Zhao  and 
(31)  Virkar  [22]. 


2.6.2.  Rational  efficiency  (ifr) 

To  give  a  realistic  indicator  of  the  system  efficiency  by  taking 
exergy  of  the  working  fluid  at  the  outlet  that  is  obviously  useful 
in  the  next  component,  the  rational  efficiency  or  rational  exergetic 
efficiency  is  applied.  The  rational  efficiency  is  one  form  of  exergetic 
efficiencies  initially  defined  by  Kotas  [30]  as  a  ratio  of  the  desired 
exergy  output  to  the  total  exergy  input.  The  rational  efficiency  can 
be  applied  to  SOFC  as  follows: 

xh  =  _ _  (32) 

^  (EXfuel  +  Exair)in  ~  i^fuel  +  £*azr)out 

2.7.  Numerical  result  validation 

To  reveal  performance  limits  and  defects,  numerical  modeling 
has  become  a  valuable  tool  for  design  and  analysis.  For  a  reliable  dis¬ 
cussion,  it  is  essential  to  validate  the  numerical  results.  Because  the 
detailed  information  on  experimental  results  for  direct-internal- 
reforming  SOFCs  can  rarely  be  found  in  literature,  in  this  study, 
the  model  verification  was  achieved  by  comparing  the  numerical 
results  with  the  actual  case  studies  and  simulation  results  by  Flo 
et  al.  [16]. 

The  electrochemical  model  used  in  this  study  was  validated  by 
comparing  the  numerical  results  with  the  experimental  data  of 
Zhao  and  Virkar  [22].  The  operating  parameters  and  cell  geometries 
are  summarized  in  Table  6.  As  shown  in  Fig.  3,  the  predicted  results 
fairly  agree  with  the  experimental  counterparts.  In  particular,  a 
good  agreement  is  achieved  when  the  operating  temperature  is  rel¬ 
atively  high  at  1073  K.  It  corresponds  to  the  temperature  range  of 
the  numerical  simulations  discussed  in  the  following  chapter;  con¬ 
sequently,  the  numerical  model  reliability  can  be  established.  The 
thermal  model  was  also  verified  comparing  the  predicted  results 
with  the  CFD-based  simulation  results  of  Flo  et  al.  for  the  planar 
anode-supported  SOFC  with  direct  CFI4  reformation  [16].  The  chan¬ 
nel  geometries  and  operating  conditions  are  illustrated  in  Table  7. 
In  spite  of  the  rather  simple  modeling  of  the  present  study  com¬ 
pared  to  the  CFD-based  2-D  model  by  Ho  et  al.,  the  performance 
predicted  by  the  present  1-D  model  agrees  reasonably  well  with 
the  results  of  Ho  et  al.  [  1 6]  as  summarized  in  Table  8.  Present  results 
show  slightly  higher  temperature  deviations  compared  to  those 
predicted  by  Ho  et  al.  for  co-flow  operation  and  lower  temperature 
deviations  for  counter-flow  operation.  This  discrepancy  between 
the  two  simulations  is  likely  ascribed  to  the  different  model  for  the 
activation  overpotential  as  well  as  the  material  properties  taken 
from  different  sources. 


3.  Results  and  discussion 

3.1.  Fundamental  characteristics  ofDIR-SOFC 

The  model  described  in  the  previous  section  is  applied  to  a  pla¬ 
nar  DIR-SOFC  to  find  its  fundamental  characteristics  and  to  confirm 
the  model’s  capability.  Considering  the  periodic  structure,  it  is  mod¬ 
eled  with  one  channel  region  of  a  single  cell  as  illustrated  in  Fig.  1, 
The  cell  is  a  typical  anode-supported  structure  with  material  prop¬ 
erties  described  in  Tables  1-5.  The  same  geometry  of  one  channel 
region  is  used  for  both  co-  and  counter-flow  cases.  For  the  sake  of 
consistency  and  simplicity,  the  inlet  gas  composition  is  chosen  as 
a  typical  CH4  syngas  composition  with  steam  to  carbon  ratio  of  2 
[31  ].  The  fuel  flows  from  left  to  right  in  the  following  figures  while 
the  air  stream  direction  is  varied  corresponding  to  the  co-flow  and 
counter-flow  configurations.  The  fixed  operating  parameters  are 
summarized  in  Table  9.  Particular  attention  is  paid  to  temperature 
profiles  of  the  PEN  layer,  since  the  cell  stack  temperature  gradi¬ 
ent  and  maximum  local  cell  temperature  are  the  most  important 
constraints  due  to  material  limitations. 

3.1.1.  Co-flow  operation 

Strong  endothermic  reaction  of  methane  reforming  proceeds 
near  the  inlet.  It  results  in  a  rapid  change  of  the  fuel  gas  composi¬ 
tions  and  temperature  dip  as  can  be  seen  in  Figs.  4  and  5.  The  steep 


Table  7 

Operating  parameters  and  SOFC  channel  geometries  used  in  model  validation  [23]. 


Parameters 

Value 

Cell  length  (mm) 

100 

Cell  width  (mm) 

1 

Air  channel  height  (mm) 

1 

Fuel  channel  height  (mm) 

0.6 

Anode  thickness  (pm) 

630 

Cathode  thickness  (pm) 

50 

Electrolyte  thickness  (pm) 

20 

Fuel  inlet  pressure  (bar) 

1 

Air  inlet  pressure  (bar) 

1 

Fuel  utilization  rate 

0.85 

Air  utilization  rate 

0.3 

Fuel  inlet  temperature  (I<) 

1073 

Air  inlet  temperature  (K) 

In  let  gas  composition 

1073 

Fuel: 

17.07%CH4,  2.40%CO,  4.91%C02) 

26.86%H2,  and  48.75%H20 

Air: 

21%02  and  79%N2 

20 
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Table  8 

Model  validation  results. 


Parameter 

Co-flow 

Counter-flow 

Ho  et  al.  [16] 

Present 

Ho  et  al.  [16] 

Present 

Voltage  (V) 

0.7 

0.69 

0.7 

0.73 

Power  density  ( W  A-2 ) 

3850 

3819 

5320 

5747 

Current  density  (iave)  (Am-2) 

Average 

5500 

5498 

7600 

7888 

Maximum 

6350 

7500 

12,400 

12,600 

Minimum 

4000 

3400 

3500 

3216 

Cell  temperature  ( TsoM )  (I<) 

Maximum 

1155 

1184 

1260 

1160 

Minimum 

1020 

1013 

1090 

1076 

Table  9 

Operating  conditions  for  SOFC  modeling. 


Parameters 

Value 

Fuel  inlet  pressure  (bar) 

1 

Air  inlet  pressure  (bar) 

1 

Fuel  utilization  rate 

0.8 

Air  utilization  rate 

0.3 

Voltage  (V) 

0.8 

Fuel  inlet  temperature  (I<) 

1073 

Air  inlet  temperature  (I<) 

Inlet  gas  composition 

1073 

Fuel: 

17.10%CH4,  2.94%CO,  4.36%C02, 

26.26%H2  and  49.34%H20 

Air: 

21%02  and  79%N2 

Fig.  4.  Fuel  channel  mole  fractions  along  the  cell  length,  co-flow  operation. 

gradients  of  the  methane  and  steam  molar  fractions  observed  in 
Fig.  4  show  their  quick  consumptions  resulting  in  the  prominent 
increase  of  the  hydrogen  molar  fraction.  Because  of  the  heat  con¬ 
sumption  associated  with  the  reforming  reaction,  a  local  minimum 


temperature  is  observed  near  the  fuel  inlet  as  shown  in  Fig.  5.  After 
methane  depleted  approximately  47  mm  from  the  fuel  inlet,  the 
electrochemical  oxidation  of  FI2  in  parallel  with  the  shift  reaction 
are  the  major  reactions  proceed  in  the  fuel  passage  and  release  reac¬ 
tion  heat.  It  raises  the  local  temperatures  of  the  gasses  and  the  cell 
resulting  in  the  maximum  cell  temperature  ( Tsoiidmax )  at  the  end 
of  the  channel  (1118  K).  The  maximum  cell  temperature  gradient 
( dTsoudldx)max  is  1.25  Kmrrr1  observed  at  the  middle  of  the  cell.  It 
is  also  noted  in  Fig.  5  that  the  fuel  flow  temperature  is  almost  same 
as  the  cell  temperature  because  of  the  small  heat  capacity  of  the 
fuel  gas  and  flow  rate.  Near  the  inlet,  the  air  temperature  is  higher 
than  the  fuel  temperature.  The  air  flow  actually  serves  as  a  heat 
source  for  the  reforming  reaction,  not  as  a  coolant,  in  this  region 
under  the  present  condition.  The  air  flow  works  as  a  coolant  only 
after  x  >  20  mm  where  the  air  temperature  is  always  less  than  the 
cell  temperature. 

Fig.  6  presents  the  distributions  of  the  open-circuit  voltage, 
polarizations,  ohmic  loss  and  local  current  density  under  the  same 
operation  condition.  The  average  current  density  was  2779  Am-2 
in  this  case.  It  shows  that  both  the  open-circuit  voltage  and  the  local 
current  density  have  maximum  values  in  the  middle  of  the  cell  but 
their  positions  do  not  coincide.  Their  distributions  are  affected  by 
the  local  temperature  and  local  gas  composition.  The  open-circuit 
voltage  is  the  highest  around  x  =  1 5  mm  where  the  value  of  the  acti¬ 
vation  overpotential  is  also  high  because  of  the  locally  reduced 
temperature.  On  the  other  hand,  although  the  activation  overpo¬ 
tential  is  the  lowest  near  the  exit,  the  open-circuit  voltage  is  also 
low  in  this  region  because  most  of  the  fuel  (hydrogen)  has  already 
been  used  up.  Consequently  the  electrochemical  reaction  is  most 
prominent  at  the  middle  of  the  cell.  The  figure  also  shows  that 
the  activation  polarization  is  the  major  loss  whereas  the  ohmic 
loss  and  the  concentration  polarization  are  relatively  low  and 
uniform. 
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Fig.  6.  Open-circuit  voltage,  polarization  terms  and  local  current  density  distribu¬ 
tions,  co-flow  operation. 
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Fig.  7.  Fuel  channel  mole  fractions  along  the  cell  length,  counter-flow  operation. 

3.1.2.  Counter-flow  operation 

Fig.  7  shows  the  distributions  of  the  mole  fractions  in  the  fuel 
stream  for  counter-flow  case.  Comparing  Fig.  7  to  Fig.  4,  it  is  noted 
that  the  reforming  reaction  proceeds  more  rapidly  in  the  counter¬ 
flow  case  and  methane  it  is  depleted  approximately  33  mm  from 
the  fuel  entrance.  This  is  caused  by  the  elevated  cell  tempera¬ 
ture  compared  to  the  co-flow  case.  The  higher  the  temperature 
is,  the  faster  the  steam  reforming  proceeds.  The  main  difference 
between  the  two  cases  is  the  air  temperature  approaching  to  the 
fuel  entrance  region,  0  <  x  <  20  mm,  where  the  reforming  reaction  is 
active.  While  the  air  flow  is  directly  supplied  to  the  active  reform¬ 
ing  region  in  the  co-flow  case,  it  goes  through  the  air  passage 
of  the  cell  in  the  counter-flow  case  accumulating  heat  generated 
in  the  cell.  Consequently,  as  shown  in  Fig.  8,  the  air  temperature 
approaching  the  active  reforming  region  becomes  much  higher  in 
the  counter-flow  case  compared  to  the  co-flow  case  even  though 
the  inlet  temperatures  are  the  same.  This  is  an  effective  heat  recov¬ 
ery  process  of  a  DIR-SOFC.  However  from  the  viewpoint  of  the 
temperature  distribution,  the  counter-flow  configuration  needs  to 
be  examined  carefully.  The  local  cell  temperature  reaches  its  max¬ 
imum  value  of  1194K  at  30  mm  from  the  fuel  entrance  and  the 
maximum  cell  temperature  gradient  (3rso^/3x)max  is  5.10  K  mm-1 
located  at  8  mm  from  the  fuel  entrance.  Comparing  to  the  co-flow 
configuration,  the  maximum  local  cell  temperature  is  increased 
76  K  and  the  maximum  cell  temperature  gradient  is  increased  8 
times.  It  should  be  noted  that  a  large  temperature  gradient  causes 
excessively  high  stress  in  the  cell  resulting  in  thermal  cracking  and 
cell  failure.  Therefore,  the  co-flow  case  is  superior  to  the  counter¬ 
flow  case  in  term  of  material  point  of  view  and  the  cell  temperature 
profile  must  be  carefully  monitored. 


Distance  from  fuel  inlet  (mm) 

Fig.  9.  Open-circuit  voltage,  polarization  terms  and  local  current  density  distribu¬ 
tions,  counter-flow  operation. 

The  open-circuit  voltage,  polarizations,  ohmic  loss  and  local  cur¬ 
rent  density  distributions  are  shown  in  Fig.  9.  The  average  current 
density  was  4394  A  m-2  in  this  case.  The  figure  shows  that  the  non¬ 
uniformity  of  the  local  current  density  distribution  is  considerably 
reinforced  in  the  counter-flow  configuration  compared  to  that  of 
the  co-flow  configuration  shown  in  Fig.  6.  The  ohmic  loss  and  con¬ 
centration  polarization  exhibit  distributions  that  track  the  local 
current  density  distribution,  whereas  raising  temperature  results 
in  the  reduction  of  activation  polarization. 

3.2.  Effects  of  current  density 

To  compare  the  performance  of  the  two  configurations,  Fig.  10 
illustrates  the  predicted  cell  efficiencies  and  power  density  as  a 
function  of  current  density  with  cell  voltage  variation  from  0.6  to 
0.8  V.  All  other  parameters  are  fixed  at  their  standard  values  in 
Table  9.  It  can  be  seen  that,  the  performance  of  the  cell  is  hindered 
with  an  increase  of  current  density.  At  1000  Am-2,  the  power  den¬ 
sity  and  efficiencies  under  the  two  configurations  are  very  close  to 
each  other.  The  rational  efficiency  is  much  higher  than  energy  effi¬ 
ciency,  the  main  exergy  losses  owing  to  internal  consumptions  and 
the  main  energy  losses  associated  with  waste  heat.  It  shows  that 
there  is  a  considerable  potential  in  SOFC  application  to  generate 
additional  electric  or  heat  power  from  the  outlet  streams.  As  cur¬ 
rent  density  increases,  the  power  density  difference  between  co- 
and  counter-flow  cases  becomes  more  pronounced.  The  energetic 
efficiencies  vary  from  60.7  to  46.9%  and  from  61.6  to  54.0%;  ratio¬ 
nal  efficiencies  vary  from  83.8  to  66.1%  and  from  84.6  to  75.3%, 
for  co-flow  and  counter-flow  cases,  respectively.  The  increase  of 
the  difference  in  the  efficiencies  between  the  co-  and  counter-flow 
cases  is  ascribed  to  the  cell  temperature  difference.  Because  the  fuel 
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Fig.  8.  Temperature  distributions,  counter-flow  operation. 


Fig.  10.  Comparison  of  efficiencies  and  power  density  versus  average  current  den¬ 
sity  for  co-  and  counter-flow  operation. 
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Fig.  11.  Cell  voltage  and  polarizations,  co-flow  operation. 


Fig.  12.  Cell  voltage  and  polarizations,  counter-flow  operation. 


and  air  utilization  factors  are  kept  constant  as  shown  in  Table  9,  the 
flow  rates  are  tuned  as  the  average  current  density  is  varied.  At  a 
high  average  current  density,  the  flow  rates  are  increased  and  the 
amount  of  heat  absorption  by  the  reforming  reaction  becomes  large 
as  well  as  the  heat  generation  in  the  cell  associated  with  the  electro¬ 
chemical  reaction  and  other  irreversible  losses.  While  the  enthalpy 
of  the  inlet  air  flow  is  immediately  supplied  to  the  reforming  reac¬ 
tion  near  the  entrance  region  in  the  co-flow  configuration,  the  air 
temperature  first  increases  in  the  counter-flow  configuration  accu¬ 
mulating  the  heat  generated  in  the  cell.  It  results  in  a  difference  of 
the  cell  temperature,  even  though  the  sum  of  heat  generation  and 
absorption  in  a  cell  is  expected  to  be  similar  in  the  two  configu¬ 
rations.  The  difference  of  the  average  cell  temperature  between 
two  configurations  becomes  larger  as  the  average  current  den¬ 
sity  is  increased.  Fig.  10  also  shows  that  the  energetic  and  rational 
efficiencies  exhibit  similar  declining  trends.  This  implies  that  the 
increases  of  the  differences  between  exergy  contents  of  inlet  and 
outlet  streams  of  the  cell  under  co-  and  counter-flow  operations  are 
proportional  to  the  increases  of  energy  in  fuel  inputs.  The  decrease 
of  cell  performance  with  increasing  the  average  current  density  is 
mainly  associated  with  the  increase  of  the  activation  polarizations. 
Fig.  10  and  the  discussion  in  the  previous  sections  lead  to  some 
considerations:  (i)  the  increase  of  losses  caused  by  the  increased 
average  current  density  results  in  the  decline  of  energetic  and  ratio¬ 
nal  efficiencies;  (ii)  co-flow  operation  is  favourable  for  operation  at 
a  low  current  density  mode  due  to  the  high  efficiency  and  smaller 
temperature  gradient. 

The  analysis  of  Fig.  1 0  shows  the  advantages  of  using  the  present 
1  -D  model  in  a  system  analysis.  Unlike  a  lump  model  that  takes  only 
the  energy  balance  into  account,  the  present  model  considers  heat 
transfer  phenomena  in  the  cell  and  can  capture  the  performance 
difference  caused  by  the  temperature  non-uniformity. 

3.3.  Effects  of  inlet  temperatures 

The  inlet  temperatures  of  the  fuel  and  air  flows  are  design 
parameters  of  a  system.  It  affects  the  temperature  distribution 
through  heat  transfer  phenomena  in  the  SOFC  and  eventually 
affects  its  performance  as  explained  in  the  previous  sections.  A 
performance  comparison  was  made  with  co-  and  counter-flow  cell 
configurations  operated  at  different  air  and  fuel  inlet  temperature 
cases  as  listed  in  Table  10.  Case  I  is  the  base  case  having  the  inlet 
fuel  and  air  temperatures  shown  in  Table  9.  In  Case  II,  only  the  inlet 
air  temperature  is  reduced  by  100K  from  the  base  case  and  kept 
at  973  K,  while  both  the  fuel  and  air  inlet  temperatures  are  set  at 
973  K  in  Case  III.  The  average  current  density  is  fixed  at  4000  A  m-2 
for  all  cases.  The  cell  terminal  voltage  and  various  losses  are  shown 
in  Figs.  11  and  12,  for  the  co-  and  counter-flow  configurations, 


respectively.  Although  the  open  circuit  voltage  increases  with 
decreasing  operating  temperature,  the  cell  terminal  voltage  shows 
the  opposite  trend  accounting  of  the  increase  of  both  ohmic 
and  activation  losses.  In  particular,  the  increase  of  the  activation 
polarizations  is  significant  among  the  three  types  of  losses.  The 
concentration  polarizations  can  be  considered  as  minor  importance 
and  are  not  significantly  affected  by  the  inlet  temperatures. 

The  cell  temperature  distributions  of  the  three  cases  for  the 
co-flow  configuration  are  shown  in  Fig.  13.  The  cell  temperature 
generally  decreases  when  the  inlet  flow  temperature  is  reduced. 
The  cooling  effect  of  the  air  temperature  is  more  prominent  com¬ 
pared  to  that  of  the  fuel  flow  mainly  because  of  its  higher  flow 
rate.  The  distribution  profile  is  generally  similar  each  other  but  the 
temperature  gradient  near  the  fuel  entrance  seems  to  be  larger 
in  Case  II  than  other  cases.  Fig.  14  shows  the  temperature  distri¬ 
butions  for  the  counter-flow  configuration.  Cooling  effects  by  the 
reduced  inlet  temperature  are  obvious  in  the  figure.  As  can  been 
seen  in  Fig.  14,  the  Case  II  shows  the  most  uniform  temperature 
distribution  among  the  three  cases. 

The  performance  comparison  results  along  with  the  maximum 
local  cell  temperatures  (TSOhd,max)  and  the  maximum  cell  temper¬ 
ature  gradients  (3TS0/ld/3x)max,  well  known  as  the  most  important 
operational  constraints  for  the  planar  SOFC,  are  listed  in  Table  10. 
It  shows  that  the  cell  performance  can  be  improved  by  setting  the 
inlet  temperatures  high  because  the  entire  cell  is  maintained  at 
high  temperature  as  shown  in  Figs.  11  and  12. 

From  the  view  point  of  a  safe  operation,  cell  temperatures  ( TsoUd ) 
and  cell  temperature  gradients  (3TsoM/3x)  must  be  monitored  with 
caution.  A  steep  temperature  gradient  and  high  cell  temperature 
can  cause  a  severe  adverse  effect  on  the  life  span  of  the  fuel  cell. 


Fig.  13.  Comparison  of  cell  temperature  distributions  based  on  case  studies,  co-flow 
operation. 


Table  10 

Summary  of  SOFC  performance  based  on  case  study. 
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Inlet  stream  temp. 

Co-flow 

Counter-flow 

Case 

TfueUn  (K) 

Tair,in  (K) 

Tsolid.max  (K) 

( dTsolid/dx)max 
(Kmnrr1) 

0 

f 

Tsolid,max  (K) 

( dTsoMldx)max 
(Kmm-1) 

1 1 

I 

1073 

1073 

1120 

1.28 

54.7 

76.7 

1192 

5.17 

58.6 

81.3 

II 

1073 

973 

1044 

3.24 

45.6 

64.0 

1077 

1.86 

52.0 

72.8 

III 

973 

973 

1021 

1.20 

41.7 

57.6 

1071 

2.87 

51.5 

70.5 

Table  11 

Performance  of  SOFC  with  thick  anode  (anode  thickness  = 

1000  |mm). 

Inlet  stream  temp. 

Co-flow 

Counter-flow 

Case 

Tfuelin  (I<) 

Tair,in  (K) 

Tsolid.mcvc  (K) 

( dTsolid/dx)max 

0 

Tsolid.max  (I<) 

( dTsoMldx)max 

1 1 

f 

(Kmm-1) 

(Kmm-1) 

I 

1073 

1073 

1116 

1.23 

54.2 

75.9 

1168 

3.39 

57.6 

79.3 

II 

1073 

973 

1041 

1.99 

45.3 

63.4 

1062 

1.40 

51.0 

70.8 

III 

973 

973 

1017 

1.05 

41.4 

57.0 

1052 

1.93 

49.7 

67.7 

Their  allowable  values  depend  on  materials,  cell  structure  and 
manufacturing  process.  In  this  study,  the  values  of  the  maximum 
allowable  temperature  gradient  and  the  maximum  allowable  cell 
temperature  are  set  at  1300 K  and  5  Kmm-1  [4],  respectively,  to 
consider  the  viability  of  SOFC.  The  maximum  local  cell  temper¬ 
atures  presented  in  Table  10  do  not  seem  to  pose  any  problem 
to  the  cell  while  the  maximum  temperature  gradient  exceeds  the 
allowable  criteria  in  one  case. 

In  case  II  under  co-flow  operation,  the  maximum  temperature 
gradient  of  3.24  K  mm-1  occurs  near  the  fuel  entrance,  while  the 
efficiencies  are  between  those  of  the  case  I  and  III.  While,  case 
I  under  counter-flow  operation  is  incompatible  with  the  opera¬ 
tional  constraints.  Although  case  I  under  counter-flow  operation 
has  the  highest  efficiencies,  it  is  not  a  preferable  operating  condition 
due  to  the  considerably  large  maximum  cell  temperature  gradi¬ 
ent  of  5.1 7  K  mm-1,  a  critical  operating  condition  for  the  cell.  By 
reducing  air  stream  inlet  temperature,  case  II  under  counter-flow 
operation,  the  maximum  cell  temperature  gradient  and  maxi¬ 
mum  local  cell  temperature  are  brought  down  to  1.86  Kmm-1  and 
1077  K,  respectively.  The  study  shows  that,  considering  the  bal¬ 
ance  between  the  cell  performance  and  the  operation  safety,  the 
most  favourable  operating  condition  among  the  six  cases  is  case  I 
under  co-flow  operation  with  the  energetic  and  rational  efficiencies 
of  54.7,  and  76.7%,  respectively,  at  moderate  maximum  cell  tem¬ 
perature  gradient  and  moderate  maximum  local  cell  temperature 
of  1.28  K  mm-1  and  1120K.  For  all  these  results,  the  performance 
of  DIR-SOFC  under  co-flow  operation  has  a  good  potential  to  be 
further  enhanced  by  simultaneously  increasing  inlet  fuel  and  air 


Fig.  14.  Comparison  of  cell  temperature  distributions  based  on  case  studies, 
counter-flow  operation. 


temperatures  until  maximum  cell  temperature  and  cell  tempera¬ 
ture  gradient  approaching  the  material  constraints. 

In  order  to  show  the  importance  to  consider  material  constraints 
more  clearly,  results  of  additional  simulations  of  DIR-SOFC  with  a 
thick  anode  are  presented.  In  this  simulation,  the  anode  thickness  is 
doubled  from  its  standard  size,  500-1000  p,m.  Other  geometric  and 
computational  conditions  are  unchanged  from  the  standard  cases 
discussed  above.  The  performance  of  the  SOFC  with  thick  anode  is 
summarized  in  Table  11.  A  comparison  with  Table  10  shows  that 
the  effect  of  the  anode  thickness  on  the  energetic  and  rational  effi¬ 
ciencies  is  minor.  The  efficiencies  of  the  thick  anode  cells  are  1-2% 
smaller  than  those  of  the  standard  cells.  The  drop  of  efficiency  is 
mainly  attributed  to  the  increase  of  ohmic  loss  and  concentration 
overpotential.  The  effect  of  the  anode  thickness  on  the  temperature 
field,  on  the  other  hand,  is  significant.  The  maximum  tempera¬ 
ture  and  temperature  gradient  are  reduced  in  the  thick  anode  cells 
compared  to  those  of  the  standard  cells.  This  is  caused  by  the  reduc¬ 
tion  of  the  thermal  resistance  of  the  cell  in  the  directions  parallel 
to  the  cell  surface.  In  all  calculation  listed  in  Table  11,  the  maxi¬ 
mum  temperature  and  maximum  temperature  gradient  are  lower 
than  the  allowable  limits,  1300  K  and  5  Kmm-1.  Increasing  anode 
thickness  allows  counter-flow  cell  configuration  to  operate  within 
a  safe  operational  condition.  In  Table  11,  considering  the  balance 
between  the  cell  performance  and  the  operation  safety,  the  most 
favourable  operating  condition  among  the  six  cases  is  case  I  under 
counter-flow  operation. 

The  above  discussion  based  on  Tables  10  and  11  show  that  a 
preferred  flow  configuration  can  be  changed  depending  on  the  cell 
geometry  and  operation  conditions,  if  the  material  constraints  are 
considered.  This  is  a  result  that  can  never  be  obtained  if  the  discus¬ 
sion  is  based  only  on  the  energy  efficiency.  The  energy  efficiency 
of  a  counter-flow  case  shown  in  Tables  10  and  1 1  is  always  higher 
than  its  counterpart  co-flow  case. 

4.  Conclusions 

A  1-D  numerical  model  for  a  planar  anode-supported  DIR- 
SOFC  with  co-  and  counter-flow  configurations  was  developed  and 
validated.  The  calculations  were  carried  out  varying  the  average 
current  density  (1000-8000  Am-2)  and  the  gas  inlet  temperatures 
(973-1073  K)  while  keeping  the  fuel  and  air  utilization  factors  con¬ 
stant  at  0.8  and  0.3,  respectively.  The  two  configurations  were 
evaluated  through  energy  and  exergy  concepts  with  a  consider¬ 
ation  for  the  material  constraints.  Careful  attention  is  paid  to  the 
maximum  local  temperature  and  the  maximum  temperature  gra¬ 
dient  of  the  cell. 
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1.  The  significant  difference  is  observed  between  the  rational  effi¬ 
ciency  and  the  energetic  efficiency  in  both  configurations.  It 
shows  potential  for  additional  power  generation  utilizing  the 
exergy  in  outlet  streams. 

2.  As  the  average  current  density  is  increased,  the  efficiencies  natu¬ 
rally  decrease  but  the  tendency  is  more  prominent  in  the  co-flow 
configuration.  The  difference  of  the  efficiencies  between  the  co- 
and  counter-flow  configurations  is  very  small  at  low  current 
density  but  becomes  more  pronounced  at  high  current  density. 
The  dependency  of  the  energy  and  rational  efficiencies  on  the 
flow  configurations  is  successfully  captured  by  applying  the  1-D 
simulation. 

3.  As  a  result  of  the  combined  effects  of  heat  generation,  heat 
absorption  and  heat  transfer  in  the  cell,  non-uniform  temper¬ 
ature  distribution  is  formed.  The  flow  configuration  affects  the 
convective  heat  transport  and  plays  a  crucial  role  in  supply¬ 
ing  heat  to  the  reforming  reaction  near  the  fuel  entrance.  The 
counter-flow  configuration  generally  achieves  high  efficiencies 
but  unfavourable  in  terms  of  strong  temperature  gradient.  A  pre¬ 
ferred  flow  configuration  can  be  changed  depending  on  the  cell 
geometry  and  operation  conditions,  if  the  material  constraints 
are  considered. 

4.  To  evaluate  the  performance  of  SOFC,  the  material  constraints 
need  to  be  considered  as  well  as  the  energy  and  rational  efficien¬ 
cies.  The  developed  1-D  model  based  on  the  energy  and  exergy 
concepts  can  capture  the  temperature  distribution  affected  by 
the  control  parameters  of  the  system,  such  as  the  gas  inlet  tem¬ 
peratures,  gas  flow  rates  and  the  average  current  density.  It  is 
a  useful  tool  for  system  analysis  work  to  improve  the  system 
design  and  reliability. 
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